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Abstract 

Mesoscopic multi-terminal Josephson junctions are novel devices that provide weak 
coupling between several bulk superconductors through a common normal layer. 
Because of the nonlocal coupling of the superconducting banks, a current flow be- 
tween two of the terminals can induce a phase difference and/or current flow in 
the other terminals. This "phase dragging" effect is used in designing a new type 
of superconducting phase qubit, the basic element of a quantum computer. Time- 
reversal symmetry breaking can be achieved by inserting a 7r-phase shifter into the 
flux loop. Logical operations are done by applying currents. This removes the neces- 
sity for local external magnetic fields to achieve bistability or controllable operations. 
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Although time-domain coherent oscillations have been observed in superconducting 
charge qubits [1], the short decoherence time t v , due to the fluctuations of the 
background charges, prevents these qubits from being a good candidate for large- 
scale quantum computing. Phase qubits, on the other hand, can couple weakly 
to the background charges and therefore potentially have larger r«. To achieve a 
reasonably long r^, it is necessary to have a "quiet" [2] phase qubit — with small 
magnetic coupling to the environment, or equivalently, small inductance. A usual 
rf-SQUID can show bistability only when the inductance L of the ring exceeds 
2tt&q/I c [3], and therefore cannot be quiet. Here, I c is the Josephson critical current 
of the junction and <I>o = h/2e the flux quantum. To overcome this problem, three 
Josephson junctions have been included in a superconducting ring [4]. One of the 
three Josephson phases is fixed by the other two and the external flux, which leaves 
the SQUID with two degrees of freedom, making bistability possible even when 
L = 0. 

A four-terminal junction is described by three phase variables (see below). Connect- 
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Fig. 1. Mesoscopic four-terminal junction (a) and four-terminal SQUID (b). The 
7T-j unction is included in the ring to attain bistability. 

ing two of the terminals by a superconducting ring will fix one of the phases to the 
external flux (neglecting the inductance of the loop, L). The resulting four-terminal 
SQUID [5] will have two degrees of freedom and can exhibit bistability at small L. 
Bistability of a four-terminal SQUID made from microbridges has been observed 
experimentally [6]. As we shall see, with a mesoscopic 4-terminal junction (Fig. la) 
it is possible to have bistability even at L = 0, due to the phase-dragging effect 



A mesoscopic 4-terminal junction is shown in Fig. la. The four bulk superconductors 
are connected to each other via a 2-dimensional electron gas (2DEG) region. The 
phase of the order parameter in the i-th terminal is denoted by fa. When the 
dimensions of the 2DEG region are smaller than the superconducting coherence 
length in the banks, the total current Ii flowing into the i-th terminal depends on 
the superconducting phases fa in all the banks through [7] 



where jij are Josephson coupling constants [9], faj = fa— fa, and Ao is the supercon- 
ducting gap. We study the system at temperatures close to T = 0, 
where decoherence due to the environment is minimal. In this limit, the Joseph- 
son energy associated with the four-terminal junction is given by 



Here, Eq = hlo/e and 1$ = Tt^\2&-o/e are the Josephson energy and critical current 
for the subjunction 1-2 at T = 0, respectively. 



[7,8], . 
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Fig. 2. Contour plot of the free energy at T = and I = 0. (a) Junction with 
e = 0. (b) Junction with e = 0.03. Other parameters are 7 = 0.1, 8\ = 82 = 0.05. 

A mesoscopic 4-terminal SQUID is constructed from the 4-terminal junction by 
connecting two of the terminals via a superconducting ring (Fig. lb, ignore the ir- 
junction for the moment). We label the terminals in such a way that subjunction 
1-2 forms the bias circuit carrying current I, and subjunction 3-4 makes the flux 
loop with current J and flux <3? threading the ring. We introduce new variables by 
01,2 = + x)/2 and ^3,4 = (±0 — x)/2, implicitly setting Yl&i = 0, which is 
allowed because the overall phase is arbitrary. The phase differences 9 and (ft are 
between terminals 1-2 and 3-4, respectively. On the other hand, x is t ne overall 
phase difference between the ring and the bias circuit. It is also useful to define the 
new dimensionless parameters 



7 = (713 + 723 

e = (713 + 723 

81 = (713 - 723 

82 = (713 - 723 
« = 734/712 • 



+ 714 + 724)/7l2 

- 714 - 724)/7l2 

- 714 + 724)/7l2 
+ 714 - 724)/7l2 



In general our system has a 3D phase space ((ft, 9, %)• However, we are interested in 
the regime where k <C 7 <C 1 and L — > 0, so that the self-generated flux by the ring 
be very small (oc "floL <C 3>o)- Therefore, (j) is practically fixed by the external field 
and/or by a 7r-phase shifter inserted into the ring (see below), and we can study the 
system in the 2D phase space of (9,x)- 

Applying an external flux <3? e = $0 /2 to the superconducting ring makes the system 
bistable (as in the rf-SQUID or 3-junction cases), meaning that the free energy of 
the system has two local minima, corresponding to opposite directions of current 
in the ring. Note that the external flux is not used to manipulate the flux (qubit) 
state. It therefore can be fixed to $o/2 for all qubits. This opens the possibility of 
replacingthe external fluxes by a -zr-phase shifter [10] in each qubit's superconducting 
ring, as shown in Fig. lb. The net effect is the same but this has the advantage 
that the 7r-phase shifter does not bring in extra coupling to the electromagnetic 
environment. In the regimeof interest (L — > 0), (ft = ir and the free energy of the 
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Fig. 3. (a) Contour plot of the free energy for the system of Fig. 2b at X = 0.05. (b) 
Solid line is the energy bias as a function of the transport current, normalized to the 
barrier height, for the system of Fig. 2b. Dashed line is the linear approximation. 

system is given by 

U = -X6 + £ J (9, X ,<f> = ir), (4) 



where I = I/Iq. Contour plots of this free energy at two different sets of parame- 
ters are given in Figs. 2 and 3a. Tunneling between the potential wells is enabled 
by charging effects, with the electrostatic capacitance of the system defining the 
effective "mass" in the kinetic energy term (see e.g. [11], sec. 2. 2. 2). 

When X = 0, the two minima of IA have equal energy. Contour plots of the free energy 
of the system at X = are shown in Fig. 2. As is clear from the figure, with the 
parameters chosen, the minima are located very close to 9 = 0. In extended phase 
space, there are also other minima, near 9 = 2im (n an integer). Those minima are 
separated from the ones shown in the figure by haigh and wide potential barriers. 
Therefore, tunneling in those directions is negligible. The situation is different for 
X- When e = 0, as is the case for a system with a square 2DEG region and four 
equivalent terminals, the minima are equidistant at x = 0, ±7r, with equal barriers 
between them (Fig. 2a). Therefore the tunneling probabilities in the left and right 
directions are the same. This is undesirable for qubit application because it makes 
the system sensitive to random charges in the environment [12,13]. However, making 
e / will move two of the minima closer together, making the barrier heights 
unequal (Fig. 2b). Pairs of minima are then isolated, and one can associate a given 
pair of minima with the logical qubit states {|0), |1)}. This regime can be achieved 
easily by choosing a rectangular 2DEG region instead of a square one [14]. 

Applying a nonzero transport current X moves the minima from being centered 
around 9 = to some 9 = 9q(I). More importantly, it removes their degeneracy. 
Fig. 3a displays the contour plot for IA using the parameters of Fig. 2b, but with 
X = 0.05. As a result of the applied current, the two minima are now clearly unequal. 

The energy difference between the two minima e( I) is plotted in Fig. 3b. As is 
evident from the figure, this energy bias is linearly dependent on I for a relatively 
wide range of the transport current: —0.3 < X < 0.3. We can therefore approximate 
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it by e(X) = EqX, where Eq is given by [14] 

(7ft ~ e5 2 )[ 7 e(5f + 5 2 ) + S^ 2 + e 2 )] 

£ ° " 4( 7 2 + e 2 ) • (5) 



To study the quantum dynamics of this system, we need to know the capacitances 
between the terminals of the 4-terminal junction. In general, there exists a capaci- 
tance between any two terminals of the system and one has to find the component 
of the effective mass tensor along the direction of tunneling in the same way as in 
Ref. [12]. However, as is clear from Fig. 2b, the difference in 9 (and also <p) from one 
minimum to another is very small compared to that in %. The tunneling is therefore 
effectively in the x direction. Using a simplified ID model we find the tunneling 
matrix element at I = to be A ~ Hu e-^ Ub/Ec , where w = [8(t 2 + e 2 )] 1 / 4 /^^ 
is the plasma frequency at the minima, E c = e 2 /2C e fr is the charging energy, C c s 
is the effective capacitance in the direction of tunneling, and 

tt- - <>-W* _ (6) 



7+|e| + \/2(7 2 + £ 2 ) 



is the barrier height between the two nearest minima [14]. 

As mentioned above, an X / lifts the degeneracy between the lowest-energy states 
and therefore stops the coherent tunneling. This energy difference induces a relative 
phase between the logical states. Therefore, control over the transport current suf- 
fices to manipulate the effective one-qubit Hamiltonian H e ^ = A(I)a x + Eq£(I)g z . 
Entangling operations between two qubits are possible through voltage-controlled 
couplings provided by additional 2DEGs [14]. Combining the two regimes of zero 
and nonzero X and using 2-qubit coupling, it is possible to perform any quantum 
gate operations [15]. 

Most of the arguments about decoherence discussed in Refs. [3,16] carry over to 
this system. There are however also two sources of decoherence different from those 
discussed in the references. The first is decoherence caused by fluctuations of the 
transport current X. This can be reduced by increasing the internal resistance of 
the current source and working at low temperatures [3,14]. Moreover, the current 
carried by the quasiparticles through the normal region can also cause decoherence. 
As shown in [17], the quasiparticle (shunt) resistance is R qp oc T cosh 2 (E a/ 2k bT), 
where ±Ea are the energies of the Andreev bound states inside the normal region. 
To achieve (exponentially) large _R qp and therefore long , it is necessary to work at 
temperatures far below Ea- In systems with a large 2DEG region, the energy scale 
Ea is inversely proportional to the dimensions of the normal region and can be much 
smaller than the gap Ao (which determines the energy scale in tunnel junctions). 
For short junctions on the other hand (which is the case here), Ea ~ Ao cos(A0/2) 
and can be large if the phase difference Acj) is not close to tt/2. A phase-dependent 
conductance in agreement with the above picture has been observed [18]. 

In the limit studied in this paper, the time scale of the dynamics is set by the 
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Josephson and charging energies, as well as by the coupling coefficients (3). For a 
junction size of 100 nm, we estimate Iq ~ 10~ 7 A [19] and C c g ~ 10~ 13 F. Taking 
7 = 0.1, S± = 82 = 0.05 and e = 0.04, we obtain A ~ 0.1GHz while tunneling 
through the barrier separating the pairs of minima is 10~ 3 smaller. The dynamics 
is thus effectively restricted to one pair of minima in phase space. Moreover, from 
(5) we obtain EqEq ~ 0.01GHz. Using the latter result, we estimate that up to 10 5 
operations can be performed within the decoherence time due to fluctuations of the 
transport current [14]. 
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